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Abstract 

The kinetic theory of Fermi liquid is applied to finite nuclei. The nuclear 
collective dynamics is treated in terms of the observable variables: particle 
density, current density, pressure etc. The relaxation processes and the de- 
velopment of instabilities in nuclear Fermi-liquid drop are strongly influenced 
by the Fermi-surface distortion effects. 
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I. INTRODUCTION. NUCLEAR FERMI-LIQUID DROP MODEL 

Dynamics and dissipative properties of nuclear Fermi liquid depend in many aspects on 
the dynamic distortion of the Fermi surface in the momentum space. It is well-known that 
taking into account this distortion allows the description of a new class of phenomena, most 
famous of which are giant multipole resonances. Furthermore, scattering of particles from 
the distorted Fermi surface leads to relaxation of collective motion and gives rise to fluid 
viscosity JT|. We point out that the development of instability in nuclear processes like 
binary fission or multifragmentation in Hi-reactions also depends on dynamic effects such 
as the dynamic Fermi-surface distortion. 

A convenient way to introduce the Fermi-liquid effects into the nuclear many body prob- 
lem is to start from the kinetic equation for the phase space distribution function f(f,p,t). 
The kinetic equation can then be transformed identically to a set (infinite) of equations for 
the moments of f(r,p, t) in p-space, namely, local single-particle density p, the velocity field 
u, the pressure tensor 7r a/ g, etc., see |l|]. In case of small variations of the particle density, 5p, 
the first order moment of the kinetic equation has the form of Euler-Navier-Stokes equation 
and is given by [§] 



The internal energy density £ in Eq. ([[]) contains both kinetic E^in and potential £ pot energy 
densities: £ = £ki n + £ po t- The pressure tensor n ai/ depends on the Fermi-surface distortion 
effect. In general case, tensor 7i aiJ also includes the viscosity tensor which is derived by the 
collision integral. 

Eq. (0) is not closed because it contains the pressure tensor 7r a p given by the second 
order moment of the distribution function f(f,p,t). We will follow the nuclear fluid dy- 
namic approach |3|-|5| and take into account dynamic Fermi-surface distortions up to the 
multipolarity of Z = 2. Assuming a periodic in time eigenvibrations with the eigenfre- 
quency uj = u + iT/2h, where Uq and Y are real, and a separable form of the velocity field, 




(1) 
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u(f,t) = (3(t)v(f), with (3{t) ~ e luJt , Eq. ([[]) is reduced to the equation of motion for the 
macroscopic variable (3{t) with the following secular equation 

-uj 2 B + (C {ld) +C'{u))-iu 1 {u) = Q. (2) 

Here, B and are the mass coefficient and the stiffness coefficient in the traditional 

liquid drop model (LDM) || respectively. The additional contribution from C'(oj) to the 
stiffness coefficient and the dissipative term 7(0;) depend on the relaxation time r and are 
given by 

nl( \ f p Oo^) 2 fdv a dv v 2 dv x \dv a 

and 

where P eq is the equilibrium pressure of the Fermi gas. The additional contribution from 
C'ioS) to the stiffness coefficient in Eq. (0) is absent in the LDM, i.e. in the liquid drop 
limit lot — > 0, and represents the influence of the dynamic Fermi-surface distortion on the 
conservative forces in the Fermi system. Finally, the dissipative term 7(0?) appears due to 
the interparticle scattering from the distorted Fermi surface. 

In general, both, C'(u) and depend implicitly on the temperature, T, via the 

dependence of the relaxation time r and of P eq on T. In cold nuclei, in the zero-sound 
limit lot — > 00, the main contribution to the stiffness coefficient in Eq. (|2]) is due to the 
Fermi-surface distortion effect given by C'{u). In Fig. 1, this effect is shown in a transparent 
way for isoscalar quadrupole excitations. As it is seen from this figure, the Fermi-surface 
distortion effect leads to a significant upward shift of the energy of vibrational states to 
locate it in the region of the quadrupole giant resonance (solid line). 

II. GIANT MONOPOLE RESONANCE AND NUCLEAR INCOMPRESSIBILITY 

We will discuss the Fermi-surface distortion effect in more detail for the case of isoscalar 
giant monopole resonances (ISGMR). This particular case is important for understanding 
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the nature of nuclear incompressibility. We will consider below the model for a Fermi-liquid 
drop having a sharp surface of the equilibrium radius Ro and the bulk density p . 
The particle density variation 5p is then given by 



6p(r, t) = r)(r, t) p #(#0 - r) + q(t) p RqS(R - r) 



(5) 



where the unknown functions 77 (r, t) and q(t) are related to each other by the particle number 
conservation and the bulk density parameter rj(r,t) is found from the equation of motion 
derived by Eq. (JJ). Namely, 



d 2 1 

m — ri = -K' V 2 7? with K' = K + K 
ot 2 9 



A'- 



(6) 



Here K is the static incompressibility 



K = R 



2 * 2 E/A 



5R 2 



R=R 



(7) 



The additional contribution to the incompressibility K' in Eq. is due to the dynamic 
Fermi-surface distortion effect |3|,[7] . The value of depends on the Landau scattering 
amplitude F . In nuclear case, F Q ~ 0, one has w 2 K. 

An essential property of a finite liquid drop having a free surface is that the motion of 
the surface should be consistent with the motion of the liquid inside the drop. This can 
be achieved by imposing a boundary condition for the compensation of the compressional 
pressure 7r rr at the liquid surface by the pressure generated by the surface tension forces 
8P a . Finally, the eigenenergies in Eq. (Bj) are given by 



hbJr. 



h 2 K> 



Jn \9mR 2Xn) 
where x n = k n R are derived from the following boundary condition: 

Xn JO On) - (fa + f^jlM = 0. 



(8) 



(9) 



Here the coefficients f a and / M are related to the surface tension and the Fermi-surface 
distortion respectively and are given by 



In the general case of Fermi-liquid drop with f ^ ^ 0, the eigenfrequency u n given in Eq. 
(§) is renormalized due to two contributions associated with the Fermi-surface distortion: 
1) the direct change of the sound velocity, i.e. in Eq. ([§]) if' appears instead of if; 2) the 
change of the roots x n of the secular equation (Qj due to additional contribution from / M ^ 
in Eq. (P). These two effects work in opposite directions: if' increases uo n while / M in Eq. 
(|) decreases it. Both values of if' and x n in Eq. depend on the Landau scattering 
amplitude Fq. In the limit of zero- to first sound transition at Fq ^> 1, one has f ^ — > and 
the solution of Eqs. @ and (H) is close do the one obtained in LDM. However, for realistic 
nuclear forces F ~ 0, the general expressions ©-© with f^^O have to be used. 

In Fig. 2 we plotted the eigenenergy, E +, of the lowest monopole mode (n = 1) given 
by Eqs. (§) and (§) as a function of mass number A. The dashed line shows the result for 
Eq+ from Eq. (|j) where if' was replaced by Kldm — 190 Me V and boundary condition 
(|) with ffj, — f a — was used. The solid line shows the result of FLDM calculations. 
We adopted the following expression for K || K = if oo + if ff v4~ 1/3 with if^ = 220 Mel/ 
and if CT = -440 MeV and for p = 0.17 /m" 3 , e F = 40 MeV, r = 1.12 fm and a = 
1.2 MeV I fm 2 . So, in spite of strong renormalization of nuclear incompressibility due to 
the Fermi-surface distortion effect, the energy of the giant monopole resonance Eq+ in the 
Fermi-liquid drop is located close to the one in the traditional LDM M where the Fermi- 
surface distortion effects are absent. The compensation of the Fermi-surface distortion 
effect in nuclear incompressibility if' due to the consistent account of the same effect in the 
boundary condition is essential for the lowest mode (n = 1) only. The above mentioned 
compensation is rather small for the highest modes with n > 2. In particular, the eigenenergy 
of the double (overtone) giant monopole resonance E + (in contrast to E + ) has a significant 
upward shift with respect to E^+ M w 2 E + given by the traditional liquid drop model. 
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III. COLLISIONAL RELAXATION ON THE DEFORMED FERMI SURFACE 



The relaxation of nuclear collective motion toward thermal equilibrium has been de- 
scribed in great detail within the framework of the kinetic theory, taking into account the 
collision integral SSt(f,p,t) PJ5|J7|J^-|TT| . The collision integral 5St(r,p,t) depends on the 
equilibrium distribution function, f eq {r,p)- In case of a finite Fermi system, the equilibrium 
distribution function, f e q{r,p), contains both the diffusive layer and oscillations in p-space 
which are caused by particle reflections at the potential surface. Both of them depend on 
the distance r. This fact also leads to the r-dependence of the collisional relaxation time, 
r 2 (r,u), given by 

l/r 2 (r, u) = - J dpp 2 Y 20 8St/ J dpp 2 Y 20 8f. (11) 

The cu-dependence of r 2 (r, u) is due to the memory effects in the collision integral 0. 

The evaluation of the collisional relaxation time r 2 (r, to) in a finite Fermi system requires 
additional accuracy. The diffusive layer in the equilibrium distribution function in momen- 
tum space can lead to a spurious effect in non-zero gain- and loss fluxes of the probability 
in the ground state of the Fermi system. However, it can be shown, see Ref. [[II], that the 



general condition for the disappearance of both above mentioned fluxes is reached due to 
the occurrence of quantum oscillations in the equilibrium distribution function f eq {r,p) in 
momentum space. The r-dependent diffusive tail of the distribution function in momentum 
space leads to an increase of the collisional damping of the collective motion in the surface re- 
gion of a nucleus and thus to an increase of the isoscalar giant quadrupole resonance (GQR) 
width. However, this increase is strongly reduced due to the above mentioned oscillations 
of the equilibrium distribution functions appearing in the collision integral. As a result, the 
collisional width of the isoscalar GQR does not exceed 30-50% of the experimental value and 
agrees with the estimates of the width where the sharp Thomas-Fermi distribution function 



is used [11 



In the case of heated nuclei, both temperature and memory effects have to be taken into 
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account in the collision integral. For frequencies small compared to the Fermi energy, one 
finds 00: 



l/r 2 = {l + (hu /27rT) 2 ]/f 2 , (12) 

which is valid for T, hu <C ep and when sharp Thomas- Fermi distribution for f e q(f,p) is 
assumed. The magnitude of f 2 can be given in the following general form at low temperatures 
T « e F [1 

f 2 {T)/h = a {±) T~\ (13) 

where the quantity is for the isoscalar mode and is for the isovector mode: = 
9.2 MeV and = 4.6 MeV [|1(J . 



In general, there are several contributions to the relaxation time r. The main contribution 
r 2 arises from the above discussed interparticle collisions and the other from collisions of 
nucleons with the moving nuclear surface (one-body dissipation with relaxation time n). 
We will also include the contribution from particle emission with relaxation time Tf. Thus, 

l/r = l/r 1 + l/r 2 + l/r T . (14) 

Introduction of the one-body dissipation reflects the peculiarities of our consideration, 
namely, restricting Fermi-surface distortions to the multipolarity I < 2 does not allow us to 
take into consideration Landau damping (fragmentation width). 

The total width, T, of a collective state can be derived from Eq. with r 2 replaced by 
the total relaxation time r from Eq. (|T4|). Assuming sharp nuclear surface, the total width 



T is reduced to the following form, see Ref. [10 



T^2xhu — 2, (15) 

1 + X (uJqT) 

where x ~ (1/6) (vf/u\) 2 and u\ is the first sound velocity in nuclear matter. We want to 
note that the total intrinsic width given by Eq. ( |T5] ) has a bell-shaped form as a function 
of y = uqt. The width T is peaked around y = yo = x^ 1 ^ 2 and the maximum value of T 



is T max = TiujqX 1 / 2 . It is easy to see that y represents the crossing point of both curves 
^ (y) ~ 1/t (zero sound regime) and T^(y) ~ r (first sound regime). Due to this fact the 
condition 

uqt = Vo = (^o^)o = x~ 1/2 (16) 

can be used as the condition for the transition from long to short relaxation time regimes. 
The magnitude of the intrinsic width F decreases when the parameter ujqt exceeds xjq. We 
have yo = 2.28 for the isovector GDR at the realistic value of x = 0.19. Such a value 
of loqt can be reached at the temperature T = T tr ~ 4.5 MeV, see Eq. ([12]). If the 
relation (|16|) with x = 1 and r = r 2 is used as the condition for transition between different 
regimes of sound wave propagation then the transition between both regimes occurs at the 
higher temperature T tr ~ lOMeV. In Fig. 3 we show the intrinsic width of the giant 
dipole resonance (GDR) in the nucleus U2 Sn. As it is seen from Fig. 3, the expression 
(|l5l) (solid line) leads to a smoother behaviour of the total intrinsic width with increasing 
excitation energy as compared to the prediction of the zero-sound regime given by the 
following condition: T pa 2%/r = F\ + r 2 + T| with Ti = 2%/ti, r 2 = 2^/r 2 and = 2h/r^ 
(dashed line). Our calculation of the intrinsic width V for the isovector GDR in the hot 
nucleus 112 Sn confirms a saturation effect in the energy dependence of V. However, we 
observe a systematic deviation of the evaluated width with respect to experimental values. 
This deviation can be reduced by varying the parameter x in Eq. flTSp. An increase in the 
parameter x improves agreement with the experiment (see dotted line in Fig. 3). 

IV. BULK INSTABILITY 

In the vicinity of the equilibrium state of the nuclear Fermi-liquid drop the stiffness 
coefficients for particle density and surface changes are positive and the nucleus is stable 
with respect to the corresponding distortions. With decreasing bulk density or increasing 
internal excitation energy (temperature) the liquid drop reaches regions of mechanical or 
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thermodynamical instabilities with respect to the small particle density and experiences 
shape fluctuations and separation into liquid and gas phases. 

Let us consider small density fluctuations 5p(r,t) starting from Eq. (|1]). We will use a 
simple saturated equation of state £ = £[p] = £un[p] + £pot[p] with 



£kin[p] 



2 m 



3 /3tt 2 x2/ 3 . 1 (Vp) 2 

b(t) p 



(17) 

£ P ot[p] =t oP 2 + t 3 p 3 + t s (Vp) 2 + e c [p\. 

The potential energy £ po t[p\ has a Skyrme-type form including the surface term ~ (Vp) 2 
and the Coulomb energy density, ec[p\- 

We will again assume a sharp surface behaviour of the equilibrium particle density. 
Assuming also Sp ~ exp(— iut) and taking into account continuity equation and Eq. (PT|). 
the equation of motion (fj) can be reduced to the following form (we consider the isoscalar 
mode): 

- mu 2 5p = ( I K - \-^—{p Q / PQ )\ V 2 5p -2{f5 + t s p ) V 2 V 2 5p, (18) 
\y o 1 — iujt J 

where Pq = Pop 2 ?/^ m, (3 = (h 2 /8m) i] and K is the incompressibility coefficient 

K^6e F {l + F ) with F = -^(to + 3i 3 Po). (19) 

Let us consider volume instability regime, K < 0, and introduce the growth rate V = —iu 
(r is real, T > 0). Using Eq. (|18|), one obtains 



«i|V-C(rV-*; s g 4 , (20) 



where 



ul = (l/9m)K, C(T) = 4Ji?_ ^ „ s = ( 2 / m )(/3 + *. p ). (21) 

3 m 1 + Ft p 



The quantity ^(T) in Eq. ( |2"0"D appears due to the Fermi surface distortion effect. Equation 
( P0|) is valid for arbitrary relaxation times r and thus describes both rare and frequent 
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collision limits as well as intermediate cases. From it one can obtain the leading order terms 
in the different limits mentioned. 

(i) Frequent collision regime: r — ► 0. High temperatures. 
The contribution from the dynamic distortion of the Fermi surface, k v , can be neglected in 
this case and we have from Eqs. fl20"D and (|), 

r 2 = | Ml |V-r( 7 /m)g 2 -^g 4 , (22) 

where 7 = (4/3)rPo/po is the classical friction coefficient. In the limit of small friction, we 
obtain 
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rn 



r 2 « K| 2 q 2 - k s q A - -L g 2 J\ Ul \* q 2 _ Kg q ± (23) 



The amplitude of the density oscillations with a certain multipolarity L, Sp^r, t), grows 
exponentially if T > 0. The expression (|23"D determines two characteristic values of the wave 
number q, namely, q ma x and q cr u where the growth rate reaches the maximum of magnitude 
T max and T goes to zero (excluding q — 0) respectively, i.e., 

r = ^max at g = qmax < Qcrit-, and V = at q = q cr it- (24) 
Both q max and q crit are given by 
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_ „ , 1 ^ = Q _ (25) 

9 — Qmax 



Qcrit 1 1 K S K S + n7 ) . , 



For a saturated nuclear liquid, one has t < 0, t 3 > and i s > 0. Thus, the critical 
value q crit from Eq. (p5|) increases with a decrease of the bulk density p , see Eqs. (|T9|) and 
(p5|). The presence of friction decreases the critical value q cr it, i-e. reduces the instability. 
The existence of the critical wave number q crit for an unstable mode is a feature of a finite 
system. The growth rate T depends on the multipolarity L of the nuclear density distortion 
and on the position of qi in the interval q — -v- q cr i t The value of qi is determined 

from the boundary condition for the pressure on free nuclear surface. For given R Q the value 
of qi increases with L at L > 2. That means that if qi < q ma x instability increases with 
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L and the nucleus becomes more unstable with respect to internal clusterization into small 
pieces (high multipolarity or multifragmentation regimes) rather than to binary fission (low 
multipolarity regime). In contrast, the binary fission is preferable if q max < qi < q C rit- 

(ii) Rare collision regime: r — > oo. Cold expansion. 
Neglecting the contribution from the relaxation, one has from Eq. ( PD| ) 

4 P 

i->2 I |2 2 / 2 4 / eg /o/^\ 

T = ui g -K„g -K 8 g*, k v = - H -. 26 

3 m p eq 

The critical value g cr .i t and the value q max are given by 

2 l M l| 2 ~ K v 2 1 2 /rjy\ 

Q'crit ? Qmax r> Qcrit' \ ' ) 

K s Z 

Thus, the distortion of the Fermi-surface leads to a decrease of the critical value q cr i t i.e. 
the Fermi-liquid drop becomes more stable with respect to volume density fluctuations due 
to dynamic Fermi-surface distortion effects. 
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FIGURES 

FIG. 1. Energies of the lowest isoscalar 2 + excitations as evaluated in the traditional liquid 
drop model [6] (LDM), i.e., at C'{uj) = and in the Fermi-liquid drop model (FLDM) taking into 
account both contributions C^ LD ^ and C'(uj) in Eq. (2). The experimental data are from Ref. [15]. 

FIG. 2. Energies of the isoscalar monopole resonances in the liquid drop model (dashed line) 
and in the Fermi- liquid drop model (solid line). Experimental data from Ref. [15]. 

FIG. 3. The intrinsic width of the isovector giant dipole resonance in the nucleus 112 Sn as 
function of excitation energy U. The experimental data were taken from Ref. [16]. Solid line 
corresponds to the total width. Dashed line is the width in the zero sound approximation at 
T = 2h/r. Dotted line is the fit to the experimental data due to the parameter x = 0.78 in Eq. 
(15), see Ref. [5]. 
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